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X3 ' Abstract 
O 

On the basis of the recently proposed self-consistent theory of metal- insula- 
tor transition in strongly disordered systems, taking into account interaction 
effects, we study transition temperature T c suppression in disordered super- 
conductors for the wide disorder interval — from weakly disordered metal 
| up to Anderson insulator, induced by "Coulomb pseudogap" formation in the 

density of states. It is shown that for a number of systems this theory provides 

rather satisfactory fit of experimental data. 
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The problem of degradation of superconducting transition temperature under strong 
disordering is relatively old |I[ . It is closely connected with the question of superconductivity 
suppression due to disorder-induced metal-insulator transition |J. A number of microscopic 
mechanisms of T c suppression were proposed, such as the growth of Coulomb pseudopotential 
the influence of Coulomb corrections to the density of states etc. In the majority 
of papers only small corrections to T c due to these effects were analyzed. 

Recently we proposed a theory of metal- insulator transitions which generalize the 
self-consistent theory of localization |||| taking into account the effects of electron-electron 
interaction. This approach has allowed us to study the behavior of the generalized diffusion 
coefficient for the wide interval of disorder parameter both for metallic and insulating regions. 
These results were used in calculations of one-particle density of states with the account of 
interelectron interactions. These calculations demonstrate the formation and the growth 
of the "Coulomb pseudogap" in the density of states close to the Fermi level. In metallic 
region this behavior of the density of states corresponds to the usual square-root Altshuler- 
Aronov correction [TU|. As disorder parameter grows and system moves towards the metal- 



insulator transition this pseudogap deepens, while the effective region of square-root behavior 
diminishes, and at the point of the metal-insulator transition the density of states at the 
Fermi level becomes equal to zero — we obtain a kind of a " Coulomb gap" . In the insulating 
region, for the band of the finite width, we obtain the typical quadratic behavior of the 
density of states close to the Fermi level, reminiscent of the Coulomb gap due to Efros and 
Shklovskii [pj] , widening with the further growth of disorder. Such behavior of the density 
of states is in qualitative agreement with experiments on the number of disordered systems 
close to the metal-insulator transition U, from amorphous alloys |I^JT3|JT5|JT6|| to disordered 
single-crystals of metallic oxides, including high-temperature superconductors ||17|| . In this 
paper the results of these calculations of the density of states are used for the numerical 
study of " Coulomb gap" effects on the T c suppression for superconductors which are close 
to the metal-insulator transition. 

We shall analyze superconductivity within the framework of the simplest BCS-model. In 



the weak coupling approximation the linearized gap-equation takes the following form ||: 
A(0 = - / d£V(& 0^(0^7^(|r) A(0, (1) 

— oo 

where iV(f) - is the averaged on disorder density of states which includes the effects of 
electron-electron interaction, V(£, CO " is the effective pairing interaction. The only difference 
with the standard approach is in the account of non-trivial dependence of iV(£) on the 
electron energy £, close to the Fermi level Ep. 

In BCS theory we usually assume the existence of some effective electron- electron at- 
traction, which is determined by the balance of pairing attraction due to electron-phonon 
interaction and Coulomb repulsion. Thus we consider the effective pairing interaction in the 
following simple form: 

v(z,e) = v c (t,t')+v ph (t,e), (2) 

where K(£,£') = V C 9(E F - \Z\)6(E F - |£'|) and V ph {^') = -V ph 9(cu D - \£\)6{u, D - |£'|) - are 
the respectively the electron-electron and electron-phonon interactions, ud - is the Debye 
frequency. The constants V c > and V P h > correspond to repulsion and attraction which 
effectively operate on rather different intervals of energy: Ep ^> up>. 

After using this expression in Eq.(l) and some transformations using the even-odd prop- 
erties of the gap function A(£) we obtain: 

A(0 = [V ph 6(u D - - V C 6(E F - 0] / dgNffl-th^W) - 

<5 c 

- V c 9(Ep - / (0it/i(^r)A(0. (3) 



We look for the solution of this equation in the usual two-step form [IS 



A(«H A * (4) 

A c , ^ D < |f | < £ F , 

where A p h, A c - are some constants which are determined from the following system of 
homogeneous linear equations, which is obtained after the substitution of Eq.(4) into Eq.(3): 



{1 - (V ph - V c )N o (0)K(^-)}A ph + V c N o (0)[K(^) - K(^)]A C = 0, 



V c N o (0)K(^)A ph + {1 + V c N o (0)[K(^) - K(^)]A C = 0, (5) 



where N (0) - is the one-electron density of states of noninteracting electrons at the Fermi 
level and we introduce 

f 1 
K(0 = J dt'-thg) 



o 



N(2T c e) 



Na(0) 



(6) 



Equation for T c follows from the usual zero-determinant condition for this homogeneous 
system: 

f = rtl + AK(^)-K(^))}- 1 , (7) 

where \i* - is the Coulomb pseudopotential, \i = V c N (0) - is the Coulomb constant, A = 
V^JVo(O) - is the pairing constant due to electron-phonon interaction. In the clean limit, 
when the density of states at the Fermi level is constant, this reduces to the usual BCS- 
expression for T c . 

Equation (7) for T c was solved numerically for disorder parameter changing in wide 
interval both for metallic and insulating regions. The density of states was calculated taking 
into account lowest order corrections over electron-electron interaction 

where G R ^(p, £) = [£ — £ p ±ry — £^p^(p, " is the retarded (advanced) electron Green's 
function, E^f A )(p, £) - is "Fock" contribution to electron self- energy PJlOfl: 

f |q|<fco 

Here D(u) - is the generalized diffusion coefficient, which is determined from the following 
self-consistent nonlinear integral equation ||[|]]: 



4 



D(lo) _ { 1 D(u) r d 3 q 1 
D ~ vriVo(O) D J (2tt) 3 -iu + D(u)q 2 + 

|q|<fco 

tt jtPo J r _fq_ q^_ 

3irirN (0)J J (2n)^-i(n + uJ) + D(n + uJ)q 2 )(-in + D(n)q 2 ) 2, { } 
u |q|<feo 

where Do = Epj 3m7 - is the usual Drude diffusion coefficient, 7 = l/2r - Born scattering 
rate, r - mean free time, ko = min{pF,l~ 1 } - cut-off in the momentum space, pf - Fermi 
momentum, / - mean free path. The data on static conductivity used below were also 
obtained by numerical solution of Eq.(lO) |||7[] . 

In Fig.l we show the behavior of the density of states close to the Fermi level which 
demonstrates the evolution of the "Coulomb pseudogap" as disorder grows. This behavior 
obviously leads to superconducting T c suppression as the system moves towards the metal- 
insulator transition. 

Fig. 2 demonstrates T c suppression as disorder parameter (pj?/) -1 grows for different val- 
ues of the Coulomb constant /1 and fixed value of pairing constant A. For large \i and growing 
disorder {pfI) 1 the value of T c drops rather fast and becomes zero in metallic region far 
enough from metal-insulator transition. For smaller values of /1 this drop of T c with growing 
disorder (p^O -1 becomes slower and for small // and large enough A (dashed curves) we get 
the possibility of superconductivity persisting even in the insulating region 0. This possi- 
bility is clearly demonstrated at the insert in Fig.2, where we show the T c dependence on the 
static conductivity of the system a for the appropriate values of A and /1. For large values 
of \i as conductivity a drops T c also drops and superconductivity is completely suppressed 
rather far from the metal-insulator transition. For small values of \x this drop of T c with a 
becomes slower and for sufficiently large values of A (dashed curves) T c remains finite even 
in the case of a — > 0. 

Fig. 3 demonstrates T c degradation with the growth of disorder parameter (pfI)' 1 for 
different values of the pairing constant A for the fixed value of the Coulomb constant fi. 
For small A and disorder parameter (ppl)^ 1 growing the value of T c drops rather fast and 
becomes zero in the metallic state far from the metal-insulator transition. As A grows this 
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drop of T c becomes slower and for large enough values of A superconductivity is completely 
suppressed only somewhere in the insulating region. At the insert in Fig. 3 we show the 
dependence of the Coulomb pseudopotential /x* on the disorder parameter (p^Z) -1 for the 
appropriate values of A and \x demonstrating rather insignificant growth of \x* with disorder 
(PfO -1 close to the point where superconductivity is completely suppressed. Apparently this 
behavior is natural enough because we neglect all the processes renormalizing the matrix 
element of Coulomb interaction in Eq.(2) due Anderson localization and electron-electron 
interactions which can actually lead to rather important growth of Coulomb pseudopotential 
close to the metal- insulator transition 

This kind of behavior of T c on the static conductivity o and on disorder parameter was 
observed experimentally in a number of disordered systems, which remain superconducting 
close to the disorder induced metal- insulator transition [l|], 0, |L2| - [[17]], Jl9| - pi] . Our 



results agree rather well with experiments on amorphous alloys of InO x [14], Nb x Si\ 



rang, Au x si^ x ra-pr 



The authors of Ref. [Q had presented the results of the measurements of disorder pa- 
rameter (p_fO _1 f° r the amorphous alloy of InO x , as well as the data for T c and static 
conductivity close to the metal- insulator transition. According to our work |5],[7|] the static 
conductivity close to the metal-insulator transition can be expressed as: 

<T = <T [(p F l)WM-l], (11) 

where o"o - is some characteristic scale of conductivity close to the metal-insulator transition, 
W c (n) - the value of disorder parameter (depending on the Coulomb constant) corresponding 
to the point of metal- insulator transition. Approximating the experimental data for InO x 
by Eq.(ll) allows us to estimate the characteristic conductivity scale o"o and also, from the 
value of W c , the Coulomb constant \i. Satisfactory correlation (Cf. the insert in Fig. 3) are 
obtained for the following values: ctq — 324.95 Q^ 1 ■ cm^ 1 , and W c ~ 0.606 giving fi ~ 1.0. 

Fig. 3 demonstrates the comparison of our results with experimental data on T c depen- 
dence on static conductivity a for the amorphous InO x using the value of T co = 3.41 K, 



u D = 112 K and E F = 9.98 • 10 4 K, [uj d / E F ] ~ 1.1 • 10~ 3 - for pure In and the given above 
values of <To and fi, which allows to estimate the pairing constant A. Satisfactory agreement 
is obtained for A ~ 0.45. Dashed curves correspond to the values of A ~ 0.4 and 0.5. 

Let us discuss now the results for T c and static conductivity dependence on the Si content 
for amorphous alloys of Nb x Sii_ x | ]T5] , P^[ and Au x Sii^ x [|T^-|2"1"|] close to the metal-insulator 



transition. Assuming that the disorder parameter in this case is just proportional to Si 
concentration, so that (pW) 1 ~ 1 — x, we can express Eq.(ll) for static conductivity in the 
following form: 

° = °* X —^, (12) 
1 — x 

where x c - the critical concentration of Nb or Au at the point of metal-insulator transition. 

Approximating the experimental data for conductivity in Nb x Sii_ x Au x Sii_ x by Eq.(12) 

allows us to estimate <j$ and critical concentration x c . Satisfactory correlation (Cf. inserts 

in Fig. 5 and Fig.6) is obtained for: 

Nb x Si 1 . x : (x - 1963.9 Q- 1 ■ cm- 1 , x c ~ 0.115; 

Au x Si^ x : a ~ 2782.13 ft" 1 • cm" 1 , x c ~ 0.14. 

Fig. 5 and Fig.6 present the comparison our results with the experimental data for T c - 
dependence on conductivity a for amorphous Nb x Si\- x and Au x Sii- x , using for the pure Nb: 
T co = 9.26 K, uj d = 276 K and E F = 6.18 • 10 4 K, [u D /E F ] ~ 3.0 ■ 10~ 3 ; while for Au x Sh_ x 
we assume T co = T cmax ~ 0.86 K, uj d = 170 K and E F = 6.42 • 10 4 K, [uj d /E f ] ~ 0.9 • 10~ 3 
with the mentioned above values of <Jq, which allows us to estimate the pairing constant 
A. Assuming for these systems the Coulomb constant /i ~ 1, the satisfactory agreement is 
obtained in case of Nb x Si\- x for A ~ 0.54 and for Au x Si\- x with A ~ 0.62. 

Surely the results presented above are essentially based upon the simplest BCS-model and 
are probably oversimplified. More rigorous approach to calculations of T c must be based upon 
Eliashberg equations and realistic models of electron-phonon interaction |18|]. Especially 
this is important in case of large enough values of A, which demonstrate the possibility of 
superconductivity persisting in the insulating phase. At the same time, in present paper 



we were not concerned with the problem of the genesis of the initial value of T c q in pure 
system, but were studying only the T c dependence on disorder. In this sense our results may 
be also qualitatively valid also in the case of strong-coupling superconductivity. We must 
also note that the more rigorous analysis is also needed taking into account disorder effects 
in the matrix element of Coulomb repulsion, which lead to the additional mechanism of T c 
degradation T c |^-f|] . In the present work we have only taken into account pseudogap effects 
in the density of states. It is possible that rather satisfactory agreement with experiments can 
signify the dominating role of pseudogap formation effects in the problem of T c degradation 
under disordering, which was claimed (on the level of small corrections) already in Ref. || . 

This work was performed with the partial support of Russian Foundation of Basic Re- 
search under the grant No. 96-02-16065, as well as under the Project IX.I of the Program 
" Statistical Physics" of the Russian Ministry of Science. 
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Figure Captions: 

Fig.l. Density of states in the case of band of finite width 2E F for ^/i = 1.0 for the 
different values of disorder parameter (j»p/) _1 : 1 - 0.1,..., 5 - 0.5 - in metallic region, 7 - 
0.7,..., 10 - 1.0 - in insulating region. Dashed curve - 6 corresponds to the point of the 
metal-insulator transition. Energy e is in the units of D kl- 

Fig. 2. T c degradation as a function of disorder parameter (pfO 1 for the fixed value of 
the pairing constant A (A = 0.5 - full curves, A = 1.0 - dashed curves) for the different values 
of Coulomb constant -^/i: 1 - 0.2, 5 - 1.0. At the insert: T c dependence upon the static 
conductivity a for the appropriate values of pairing constant A and Coulomb constant /i. 

Fig. 3. T c degradation as a function of disorder parameter (pfO -1 f° r the fixed value of 
Coulomb constant = 0.4 and different values of the pairing constant A: 1 - 0.3, 2 - 0.4, 

8 - 1.0. At the insert: dependence of Coulomb pseudopotential ji* on disorder parameter 
fo r the appropriate values of pairing constant A and Coulomb constant fi. The arrow 
shows the point of the metal-insulator transition. 

Fig. 4. T c dependence on conductivity a for the amorphous alloys of InO x . At the insert: 
approximation of conductivity dependence upon the disorder parameter (pfO 1 - 

Fig. 5. T c dependence on conductivity a for the amorphous alloys of Nb x Sii^ x . At the 
insert: approximation of conductivity dependence upon concentration of Nb. 

Fig. 6. T c dependence on conductivity a for the amorphous alloys of Au x Si\- x . At the 
insert: approximation of conductivity dependence upon concentration of Au. 



9 



REFERENCES 

[1] Belitz D., Kirkpatrik T.R. Rev. Mod. Phys., 66, 261, (1994). 

[2] Sadovskii M.V. Superconductivity: Physics, Chemistry, Technology, 8, 337, (1995); 
Physics Reports, (1996), (in press). 

[3] Anderson P.W., Muttalib K.A., Ramakrishnan T.V . Phys. Rev., B28, 117, (1983). 

[4] Bulaevskii L.N., Sadovskii M.V. J.Low-Temp.Phys., 59, 89, (1985). 

[5] Belitz D., at all. Phys. Rev., B40, 111, (1989). 

[6] Kuchinskii E.Z., Sadovskii M.V., Suvorov V.G., Erkabaev M.A. 
Zh.Eksp.Teor.Fiz.(JETP), 107, 2027, (1995). 

[7] Kuchinskii E.Z., Erkabaev M.A. Fiz.Tverd.Tela (Solid State Physics), (1997), (in press). 

[8] Vollhardt D., Wolfle P. In "Electronic Phase Transitions." 

Ed. by Hanke W., Kopaev Yu., North- Holland, Amsterdam, p.l, (1992). 

[9] Sadovskii M.V. Sov. Sci. Rev., Phys., A7, 1, (1986). 

[10] Altshuler B.L., Aronov A.G. Zh.Eksp.Teor.Fiz.(JETP), 77, 2028, (1979). 

[11] Efros A.L., Shklovskii B.I. - In "Electron-Electron Interaction in Disordered Systems." 
Ed. by Pollak M., Efros A.L., North-Holland, Amsterdam, p.409, (1984). 

[12] McMillan W.L., Mochel J. Phys. Rev. Lett., 46, 556, (1981). 

[13] Imry Y., Ovadyahu Z. Phys. Rev. Lett., 49, 841, (1982). 

[14] Fiory A.T.,Hebard A.F. Phys. Rev. Lett., 52, 2057, (1984). 

[15] Hertel G., Bishop D.J., Spencer E.G., Rowell J.M., Dynes R.C. 
Phys. Rev. Lett., 50, 743, (1983). 

[16] Bishop D.J., Spencer E.G., Dynes R.C. Sol. Stat. Elect., 28, 73, (1985). 

10 



[17] Srikanth H., Rajeev K.P., Shivashankar G.V., Raychaudhuri A.K. 
Physica, C195, 87, (1992). 

[18] Vonsovsky S.V., Izyumov Yu.A., Kurmaev E.Z. Superconductivity of Transition Metals, 
Their Alloys and Compounds, (Springer Series in Solid State Sciences, Vol.27), Springer 
Verlag, 1982. 

[19] Nishida N., Yamaguchi M., Furubayashi T., Morigaki M., Ishimoto H., Ono K. 
Sol. Stat. Comm., 44, 305, (1982). 

[20] Furubayashi T., Nishida N., Yamaguchi M., Morigaki M., Ishimoto H. 
Sol. Stat. Comm., 55, 513, (1985). 

[21] Nishida N., Furubayashi T., Yamaguchi M., Morigaki M., Ishimoto H. 
Sol. Stat. Elect., 28, 81, (1985). 



11 



